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Abstract We have carried out 12 CO J=2-l and 12 CO J=3-2 observations toward the 
high-mass protostellar candidate IRAS 20188+3928. Comparing with the previous obser- 
1 vations, the 12 CO J=2-l and 12 CO J=3-2 lines both have asymmetric profiles with an 

absorption dip. The velocity of the absorption dip is ~ 1.0 km s _1 . The spectral shape 
may be caused by rotation. The velocity-integrated intensity map and position-velocity 
diagram of the 12 CO J=2-l line present the obvious bipolar component, further verifying 
this region has an outflow motion. This region also is associated with an HII region, an 
IRAS source, and H2O maser. The H2O maser has the velocity of 1.1 km s _1 . Comparing 
with the components of the outflow, we find that the H2O maser is not associated with the 
outflow. Using the LVG model, we obtained that possible averaged gas density of the 
blueshifted lobe and redshifted lobe are l.Ox 10 5 cm -3 and 2.0x 10 4 cm -3 , while kinetic 
temperature are 26.9 K and 52.9 K, respectively. Additionally, the outflow has the higher 
integrated intensity ratio (I C o .7=3-2! Ico j= 2 -i). 
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1 INTRODUCTION 

IRAS 20188+3928 is a high-mass protostellar candidate (Zhang et al . 120051 1. located in the Cygnus re- 
gion at the position {I, b) = (77 '.46, 1.76) (Little et al. [1988 ). The distance estimated to this source varies 
from 0.31 to 4 kpc (Little et al. 119881; Palla et al. I199H Molinari et al. 11996b , hence, the luminosity is 
highly uncertain. A bipolar molecular outflow was discovered in HCO + J=2-l and 13 CO J=l-0 emis- 
sion lines in the IRAS 20188+3928 region (Little et al. 1988). They suggested that the outflowing gas 
has a dense and clumpy nature in the north-SW direction, and IRAS 20188+3928 is not obviously asso- 
ciated with a well known stellar cluster and HII region. Zhang et al. (2005) observed the bipolar outflow 
in 12 CO J=2-l centered on the IRAS source and roughly in the NS direction, but their observation can 
not fully cover the emission of the bipolar outflow. In addition, Varricatt et al. (120101) not only detected 
the bipolar outflow in H2 line, but also an additional outflow in the region. Because the outflows in this 
region are complicated, farther observations are needed toward the IRAS 20188+3928 region. 

Additionally, IRAS 20188+3928 is associated with a compact molecular cloud (Anglada et al. 
I1997I ). An H2O maser emission was detected by several other investigators (Palla et al. 119911; Brand et 
al. 119941; Jenness et al. 1 1995t >. Several observers detected NH3 emission (Molinari et al. [19961 Anglada 
et al. 119971 Jijina et al. [1999), implying that the IRAS source is deeply embedded in the high density 
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gas, which may be the driving source of the bipolar outflow. Jenness et al. ( 119951 1 imaged two sources at 
450 and 800 m, agreeing with the position of two 6-cm radio emission sources detected by Molinari et 
al. (11998b . the brighter and the fainter ones located north-west of the IRAS position by 4.2 and 28.4 arc- 
sec, respectively. Near-infrared polarimetry revealed an illuminating source, about 6" north of the IRAS 
source (Yao et al. |2Q001 >. This illuminating source is better centered on the CO outflow than the IRAS 
source. Accordingly, near-IR images show a cluster of deeply embedded objects towards the centre of 
the IRAS 20188+3928 region, very close to the IRAS 20188+3928 position. The agreement of these 
positions indicates that the outflow source is located in the cluster. Hence, there are multiple YSOs in 
the region, the driving source of the bipolar outflow have not been well identified. 

In this paper, we have carried out the 12 CO J=2-l and 12 CO J=3-2 observations towards the high- 
mass protostellar candidate IRAS 20188+3928. Clearly north-SW outflow and possible rotation may be 
associated with IRAS 20188+3928. 

2 OBSERVATIONS 

The mapping observations of IRAS 20188+3928 were performed in 12 CO J=2-l and 12 CO J=3-2 lines 
using the KOSMA 3m telescope at Gornergrat, Switzerland, in November 2009. The half-power beam 
widths of the telescope at the observing frequencies of 230.538 GHz and 345.789 GHz are 130" and 
80", respectively. The pointing and tracking accuracy was better than 20". The medium and variable 
resolution acousto optical spectrometers have 1501 and 1601 channels, with total bandwidths of 248 
MHz and 544 MHz. The channel widths of 165 and 340 kHz correspond to velocity resolutions of 
0.21 and 0.29 km s _1 , respectively. The beam efficiency B c h is 0.68 at both 230 GHz and 0.72 at 345 
GHz. The forward efficiency F e g is 0.93 for all frequencies. Mapping observations were centered at 
RA(J2000)=20 h 20 m 39.30 s , DEC(J2000)=39°37'51.90" using the on-the-fly mode, the total mapping 
area is 11' x 11' with a 1' x 1' grid. The correction for the line intensities to the main beam temper- 
ature scale was made using the formula T m b = (i^ff /-Bcff x ?a). The data were reduced using the 
GILDAS/CLASS Q package. 

The 1.4GHz radio continuum emission data were obtained from the NRAO VLA Sky Survey 
(NVSS; Condon et al. 1998). 

3 RESULTS AND DISCUSSION 
3.1 CO molecular spectra 

Fig. 1 shows the spectra of 12 CO J=2-l, 12 CO J=3-2, and 13 CO J=l-0 at the IRAS 20188+3928 
position. The 13 CO J=l-0 data is obtained from the Purple Mountain Observatory (PMO) archive dat£0. 
Each spectrum display broad line wings. 12 CO J=2-l and 12 CO J=3-2 lines present asymmetric profiles 
with double peaks, while 13 CO J=l-0 line shows a single peak profile. Wu et al. (20101) only detected 
one compact core from the IRAS 20188+3928 region in the HCN J=l-0 and CS J=2-l lines, suggesting 
that there is one velocity component. The profiles of 12 CO J=2-l and 12 CO J=3-2 with double peaks 
may be a indication of absorption. The velocity of an absorption dip is ~ 1.0 km s — 1 . Rotation can 
be responsible for the spectral shape. With the single-dish data currently available, we did not find any 
kinematic evidence of rotation. The 13 CO J=l-0 line is optically thin, which can be used to determine 
the systemic velocity. A systemic velocity of ~ 2.0 km s _1 is obtained from this line. According to the 
Galactic rotation model of Fich et al. (1989) together with i? Q = 8.5 kpc and Vq = 220 km s _1 , where 
Vq is the circular rotation speed of the Galaxy, we obtain a kinematic distance of 0.32 kpc to IRAS 
20188+3928, which is used to calculate the physical parameters in this paper. The full widths (FW) of 
the 12 CO J=2-l and 12 CO J=3-2 lines both are about 20 km s" 1 from -10 km s" 1 to 10 km s^ 1 . The 
large FW appear a strong indication of outflow motion. 



1 http://www.iram.fr/IRAMFR/GILDAS/ 

2 http://www.radioast.csdb.cn 
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Fig. 1 Spectral profiles observed at the central position of IRAS 20188+3928 in the optically 
thick 12 CO J=2-l, 12 CO J=3-2 lines, and the optically thin 13 CO J=l-0 line. The dotted line 
in the spectra marks the cloud systemic velocity. 

3.2 The bipolar outflow 

3.2.1 The outflow morphology 

To determine the velocity components and morphology of the outflow, we made a position-velocity (PV) 
diagram with a cut along the north-south direction as shown in Fig. 2. The PV diagram in Fig. 2 clearly 
shows bipolar components, but we cannot see the emission from the core of IRAS 20188+3928, indi- 
cating that above cut direction represents that of the bipolar components. The blueshifted and redshifted 
components have obvious velocity gradients from —7.6 to —0.6 km s _1 and 3.2 to 9.0 km s _1 , respec- 
tively. The distributions of redshifted and blueshifted velocity components in Fig. 2 further confirm a 
bipolar outflow in this region. Using the velocity ranges of the blueshifted and redshifted components, 
we made the integrated intensity map as shown in Fig. 3. In Fig. 3, the blue and red contours represent 
the blueshifted and redshifted components of the outflow. The outflow clearly displays the north-SW 
components, which is associated with observation in HCO + J=2-l and 13 CO J=l-0 emission lines 
(Little et al. 1988). The 1.4 GHz continuum emission is superimposed on the outflow, which may come 
from an HII region (Urquhart et al. 120091 ). Both HII region and infrared source IRAS 20188+3928 are 
located the beginning position of the outflow, because of our lower spatial resolution we cannot identify 
which object is the driving source of the outflow. An H2O maser was detected in this region, which has 
the velocity of 1.1 km s _1 (Palla et al. I199H Brand et al. 119941 Jenness et al. |19951 >. Comparing the 
components of the outflow, we find that the H 2 maser is not associated with the outflow, but with the 
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Fig. 2 P-V diagram constructed from the 12 CO J=2-l transition. Contour levels are 1, 2, 
3, 4.5, 6, 6.5, 7.5, 9, 10.5, 12, 14,..., 25 K. The vertical dashed and solid lines indicate the 
beginning of the blueshifted and redshifted emission, respectively. 



core. If an accretion motion is further confirmed in this region, then the H2O maser may come from the 
disk. 

3.2.2 The outflow parameters 

In this section, we use a large velocity gradient (LVG) radiative transfer model (Goldreich et al. |1974t 
Qin et al. 2008a) to determine the gas density (rt(H 2 )) and kinetic temperature (7ki n ) of the outflow. 
The blueshifted and redshifted velocity intervals (Av) are from —7.6 to —0.6 km s" 1 and from 3.2 to 
9.0 km s _1 , respectively. In this model, assuming a uniform density and using above velocity intervals, 
n(H2) and Tk; n are obtained by fitting the line intensity of 12 CO J=2-l and the line ratio of 12 CO 
J=3-2/ 12 CO J=2-l, which are varied in a reasonable range for this region. Fig 4. shows the results 
of the LVG modeling. From the Fig. 4, various column densities are from 3.5 x 10 4 to 1.7x 10 5 cm -3 
with corresponding Tki n ranging from 22.9 to 30.9 K for the blueshifted lobe. The average values of 
possible n(H2) and Tk; n are l.Ox 10 5 cm -3 and 26.9 K, respectively. For redshifted lobe, n(H2) ranges 
from 1.6x 10 4 to 2.8x 10 4 cm~ 3 , while T kin is from 39.8 to 83.1 K. The averaged n(H 2 ) and T kin are 
2.0 x 10 4 cm -3 and 52.9 K, respectively. 

Additionally, if both lobes are approximately spherical in shape, their mass is given by 
M=n(H2)^7r_L 3 ^ig7n(H2) (Garden et al. |1991| l, where ^ g =1.36 is the mean atomic weight of the gas, 
L is the lobe diameter , and m(H 2 ) is the mass of a hydrogen molecule. We calculate the dynamic 
timescales using = r/v, where v is the maximum flow velocity relative to the cloud systemic veloc- 
ity , and r is the length of the begin-to-end flow extension for each lobe. The mass entrainment rate of the 
outflow is determined by M = M / (t^). The momentum P and energy E are calculated by P = MV 
and E = MV 2 , where V is the mean velocity of the gas relative to the cloud systemic velocity. The 
physical parameters and the calculated results of the outflow are listed in Table 1. A total mass and 
mass entrainment rate of the outflow are 3.5 M Q and 6.2 x M© yr _1 , respectively. The average 
dynamical timescale is about 6.2 x 10 4 yr. 

3.2.3 The integrated intensity ratio (Icoj=5-2/Icoj=2-i) 

The line intensity ratios based on the optically thick CO transitions can indicate the temperature varies 
at different positions (Qin et al. |2008bl , and further trace shocks (Xu et al. 1201 lal 1201 lbl [2012b . 
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Fig. 3 The velocity-integrated intensity map of 12 CO J=2-l outflow (red and blue contours) 
overlaid with the 1.4 GHz emission (black contours). The black contour levels are at 3, 6, 
9, 12, 15, 18, and 21 a (1 a is 0.002 Jy beam -1 ). The red and blue contour levels are 35, 
50,..., 95 % of the peak value. The H2O maser is shown by a filled triangle. The "*" symbols 
indicates the position of IRAS 20188+3928 and HII region. 



Table 1 The physical parameters of the outflow. 
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When the outflow comes into contact with the surrounding interstellar medium (ISM), this process 
can product the shock. In order to obtain the integrated intensity ratio of 12 CO J=3-2 and J=2-l 
(-fco ,7=3-2/^00 j=2-i), we convolved the 80" of 12 CO J=3-2 data with an effective beam size of 
Vl30 2 — 80 2 = 102". The integrated intensities were calculated for 12 CO J=2-l line in the same ve- 
locity range as for 12 CO J=3-2. Fig. 5 shows the distribution of the ratio (color scale) overlaid with the 
distribution of the 12 CO J=2-l line integrated intensity of the outflow (black solid and dashed contours). 
From the Fig. 5, we can see that the distribution of the ratios both presents shell-like morphology. The 
maximum ratio value is 1.2 in the blueshifted lobe, while the redshifted lobe has the maximum ratio 
value of 1.7. The rms level is 0.10 (la). The molecular clouds associated with HII regions have higher 
ratio and have higher gas temperature than those sources without HII regions, indicating that the high 
line ratio may be due to heating of the gas by the massive stars (Wilson et al. 1997). While the high 
ratio values in the molecular clouds interacting with SNR (Xu et al. 120 1 1 al 120 1 1 bl f20 1 21 l exceed previ- 
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Fig. 4 LVG model analysis using the 12 CO J=2-l and 12 CO J=3-2 lines. Top: The 
blueshifted components of the outflow. Bottom: The redshifted components of the outflow. 
The horizontal and vertical axes are the volume density tih 2 and kinetic temperature T^ n in 
log scale, respectively. Blue curves represent line intensity of 12 CO J=2-l, while red curves 
represent the line ratios of CO( J=3-2)/CO( J=2-l). The results of Tkin and nn 2 are taken from 
the cross points of red and blue curves. The error bars indicate the range of possible solutions. 



ous measurement of individual Galactic molecular clouds, implying that the SNR shock has driven into 
the molecular clouds. Here, the high ratio values detected in the region of the outflow may be caused 
by the outflow interacting with the surrounding ISM. In addition, the maximum ratio value of the red- 
shifted lobe is greater than that of the blueshifted lobe, we suggest that the redshifted lobe has the higher 
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Fig. 5 12 CO J=2-l intensity maps of the outflow are superimposed on the line intensity 
ratio maps (color scale). The wedge indicates the line-intensity ratio scale. Left panel: the 
blueshifted lobe; Right panel: the redshifted lobe. 



temperature, which is associated the results derived from the LVG model. In section 3.2.2, we obtain 
that the redshifted lobe has the kinetic temperature of 52.9 K, which is larger than that (26.9 K) of the 
blueshifted lobe. 



4 CONCLUSIONS 

We present submillimeter/millimeter observations of the high-mass protostellar candidate IRAS 
20188+3928 in the 12 CO J=2-l and 12 CO J=3-2 lines. The spectra profiles of the 12 CO J=2-l and 
12 CO J=3-2 both with an absorption dip, together with the 13 CO J=l-0 line with a single peak profile, 
indicate that possible rotation may be associated with IRAS 20188+3928. The 12 CO J=2-l velocity- 
integrated intensity map, position-velocity diagram, and the broad wing (full width = 20 km s _1 ) of 
the 12 CO J=2-l line, further verify that there is an outflow motion. The outflow emission is elongated 
along the north-SW direction. This region also is associated with an HII region and an IRAS source, 
but we cannot identify which object is the driving source of the outflow owing to our lower spatial 
resolution. An H2O maser is also associated with this region. Comparing with the components of the 
outflow, we find that the H 2 maser is not associated with the outflow, but may be associated with the 
core of IRAS 20188+3928. Using the LVG model, the possible kinetic temperature of the blueshifted 
lobe and redshifted lobe are 26.9 K and 52.9 K, respectively. The total gas mass, average dynamical 
timescale, and mass entrainment rate of the outflow are 3.5 M Q , 6.2 x 10 4 yr and 6.2 x 10 -5 M Q yr _1 , 
respectively. The outflow has the higher integrated intensity ratio (Ico 7=3-2/^00 J=2-i), which may 
be caused by the outflow interacting with the surrounding ISM. The maximum ratio value of the red- 
shifted lobe is greater than that of the blueshifted lobe, we suggest that the redshifted lobe has the higher 
temperature. 
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